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Abstract. In atmospheric Cerenkov technique 7— rays are 
detected against abundant background produced by hadronic 
showers. In order to improve signal to noise ratio of the ex- 
periment, it is necessary to reject a significant fraction of 
hadronic showers. The temporal and spectral differences, 
the lateral distributions and density fluctuations of Cerenkov 
photons generated by 7— ray and hadron primaries are of- 
ten used for this purpose. Here we study the differences 
in Cerenkov photon density fluctuations at the observation 
level based on Monte Carlo simulations. Various types of 
density fluctuations like the short range (or local), medium 
range fluctuations and flatness parameter are studied. The 
estimated quality factors reflect the efficiencies with which 
the hadrons can be rejected from the data. It has been found 
that we can reject around 80% of proton showers while re- 
taining about 70% of 7— ray showers in the data, based only 
on the differences in the flatness parameter. Density fluctua- 
tions particularly suited for wavefront sampling observations 
seem to be a good technique to improve the signal to noise 
ratio. 



1 Introduction 

In a typical wavefront sampling experiment, arrival time of 
Cerenkov photons and Cerenkov photon density are sam- 
pled at several locations in the Cerenkov pool generated by 
air showers initiated by 7— rays from astronomical sources, 
using distributed array of telescopes. Cosmic ray showers 
which also give rise to Cerenkov light similar to that pro- 
duced by 7— rays constitute abundant background against which 
the 7— ray signal is to be detected. Hence it is necessary 
to devise methods to reject a large fraction of cosmic ray 
background and thereby improve the signal-to-noise ratio or 
sensitivity of the experiment. Previously we have studied 
the usefulness of parameters based on timing information 
recorded in wavefront sampling experiment for gamma-hadron 
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separation (Chitnis and Bhat, 2001). Use of arrival time jit- 
ter and parameters based on shape of the Cerenkov pulse has 
been demonstrated. Here we investigate the efficacy of cer- 
tain parameters based on Cerenkov photon density distribu- 
tion for gamma-hadron separation. 

2 Local density fluctuations 

Pachmarhi Array of Cerenkov Telescopes (PACT) consists of 
25 telescopes with each telescope consisting of para-axially 
mounted parabolic mirrors of diameter 0.9 m each (see Chit- 
nis et ai, 2001 for details). Here we study the usefulness of 
local density fluctuations or LDF for gamma-hadron separa- 
tion. LDP or density jitter is defined as the ratio of RMS to 
mean of photon densities from 7 mirrors of each telescope. 
We have simulated a large number of showers initiated by 
7— rays of energy 500 GeV and protons of energy 1 TeV, 
incident vertically at the top of the atmosphere, for this pur- 
pose. Showers are simulated using CORSIKA (Heck et al., 
1998). An array of telescopes spread over an area of 400 m 
X 400 TO, much larger than PACT, is used for simulations. 
We have simulated 100 showers for each primary for each of 
the three observation altitudes, viz., sea level, 1 km above 
sea level which corresponds to altitude of PACT and for 2.2 
km above sea level. 

Figure 1 shows the variation of LDF as a function of dis- 
tance from core of the shower for showers initiated by 500 
GeV 7— rays and 1 TeV protons at various observation al- 
titudes. It can be seen that the LDF for protons is consis- 
tently higher than that for 7— ray primaries for all altitudes 
and at all core distances. This is expected due to differ- 
ences in kinematics of these two types of showers (Chitnis 
and Bhat, 1998). Also fluctuations are not very sensitive to 
core distances. Hence LDF is a likely parameter to be used 
for gamma-hadron separation. 

We use quality factor as a figure of merit of a parameter to 
distinguish between 7— ray and proton initiated showers. It 
is defined as 



2 



Bhat and Chitnis : Gamma-Hadron Separation using Density Fluctuations 



(a) 

, + ++ + + * 



+ + + ++ + + + + ^ + +++ + 



400 |- 
; 300 I- 




(°) 



0,50 1 
I 0.40 I 
; 0.30 I 

- 0.20 I 
I O-'O I 
0.00 1 


0,50 

i n.40 

] n.,"50 

^ G.?0 
: 0,10 
0.00 



^ ^ ^ J 



h + + + ++ + + + + ^. .^ooo<>o»««* 



O O O OOOOOO' 



ic) J 




(t>) J 



Cor. di.tonc, (m) 



Fig. 1. Variation of local density fluctuations (LDF) as a function 
of core distance averaged over 100 showers initiated by 500 GeV 
7— rays (indicated by diamond) and 1 TeV protons (indicated by 
+), for three different altitudes: (a) sea level, (b) 1 km from sea 
level and (c) 2.2 km. 



Fig. 2. Gamma-hadron separation based on LDF for 500 GeV 
7— rays (continuous line) and 1 TeV protons (dotted line) for obser- 
vation altitudes of (a) km, (b) 1 km and (c) 2.2 km. Dashed line 
indicates the parameter value which yields optimum quaUty factor 
(see table 1). 
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where A^^^ is the number of 7— rays accepted (i.e. below 
threshold), A'^ is the total number of 7— rays, N^^ is the 
number of protons accepted and Nj! is the total number of 
protons. Larger the quality factor, better is the background 
rejection efficiency. 

Figure 2 shows the distribution of LDF from 500 GeV 
7— rays and 1 TeV protons, for three different altitudes. Op- 
timum quality factors for each of the cases are listed in Table 
1 . It can be seen that it is possible to reject about 50% proton 
showers retaining about 85% of 7— ray showers, using LDF. 

3 Medium range density fluctuations 

PACT consists of four sectors of six telescopes each. We de- 
fine medium range or sector- wise density fluctuations (MDF) 
as the ratio of RMS to mean density, where RMS and mean 
are calculated using total photon densities at each of the six 
telescopes of a sector. Figure 3 shows the variation of medium 
range density fluctuations as a function of core distance for 
showers initiated by 500 GeV 7— rays and 1 TeV protons for 
three different altitudes of observation level. As in the case 



of LDF using the distributions of MDF for showers initiated 
by 500 GeV 7— rays and 1 TeV protons, quality factors are 
calculated at three altitudes of observation levels. Optimum 
quality factors for different altitudes are listed in Table 2 and 
distributions of MDF are shown in Figure 4. It can be seen 
that it is possible to reject about 60-70% of proton showers 
retaining about 80% of 7— ray showers, based on MDF. 



4 Medium range flatness parameter 

Lateral distributions of Cerenkov photons (variation of Cerenkov 
photon density as a function of core distance) from showers 
initiated by 7— rays show a characteristic hump at the core 
distance of about 120-140 m, depending on the observation 
altitude. This is due to the effective focusing of Cerenkov 
photons from a large range altitudes. Distributions are flat 
within the hump region and density falls rapidly beyond the 
hump. Lateral distributions from proton showers, on the other 
hand, show continuously falling density distribution as core 
distance increases (Rao and Sinha, 1988 and Chitnis and 
Bhat, 1998). Also due to the kinematical differences, the lat- 
eral distributions from 7— ray showers are smooth compared 
to proton showers. These differences in lateral distributions 
can be parameterized using flatness parameter defined as : 



Table 1. Quality factors for local density fluctuations 

Obs Threshold Quality factor Fraction Fraction 

altitude value of of of 

(km) LDF accepted accepted 

7— rays protons 

~0 0.248 1.284 ± 0.012 0.842 0.429 

1 0.257 1.215 ±0.011 0.888 0.534 

2.2 0.234 1.261 ± 0.012 0.870 0.476 



Table 2. Quality factors for medium range density fluctuations 

Obs Threshold Quality factor Fraction Fraction 

altitude value of of of 

(km) MDF accepted accepted 

7— rays protons 

~0 0.230 1.297 ±0.030 0.819 0.399 

1 0.245 1.332 ±0.031 0.830 0.388 

2.2 0.243 1.418 ±0.034 0.795 0.315 
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Fig. 3. Variation of medium range or sector-wise density fluctua- 
tions as a function of core distance averaged over 100 showers ini- 
tiated by 500 GeV 7— rays (indicated by diamond and continuous 
line) and 1 TeV protons (indicated by plus sign and dotted line), for 
three different altitudes: (a) km, (b) 1 km and (c) 2.2 km 



Fig. 4. Gamma-hadron separation based on medium range density 
fluctuations for 500 GeV 7— rays (continuous line) and 1 TeV pro- 
tons (dotted line) for observation altitudes of (a) km, (b) 1 km 
and (c) 2.2 km. Dashed line indicates the value of the parameter for 
optimum quality factor (see table 2). 
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where N : no. of telescopes triggered, pi : photon density 
measured by individual telescopes and po : average density. 

Lateral distributions from 7— ray showers are expected to 
have a smaller value of a parameter compared to proton gen- 
erated showers, on the average. Figure 5 shows variation of 
a parameter as a function of core distance for showers gen- 
erated by 500 GeV 7— rays and 1 TeV protons, for three ob- 
servation altitudes. It can be seen that, on an average, 7— ray 
showers have smaller value of a compared to proton show- 
ers. Also for both 7— rays and protons, value of a increases 
with increase in altitude of observation. At all the three alti- 
tudes 7— ray showers show larger value of a near hump re- 
gion of lateral distribution. As a result, the difference in value 
of a for 7— rays and protons near hump is reduced. Hence 
a can be useful discriminant at core distances away from the 
hump, on both the sides. 

Distributions of a parameter for telescopes within core 
distances of 100 m, for showers initiated by 500 GeV 7— rays 
and 1 TeV protons, for observation altitude of 1 km are shown 
in Figure 6. Optimum quality factors for all the three alti- 



Table 3. Quality factors for medium range flatness parameter for 
core distance < 100 m 



Obs 


Threshold 


Quality factor 


Fraction 


Fraction 


altitude 


value of 




of 


of 


(km) 


a 




accepted 


accepted 








7— rays 


protons 





0.33 


1.800 ±0.119 


0.690 


0.147 


1 


0.52 


1.556 ± 0.095 


0.738 


0.225 


2.2 


0.66 


1.662 ± 0.108 


0.677 


0.166 



tudes are listed in Table 3. It can be seen that the flatness 
parameter serves as a good discriminant for showers with 
smaller impact parameters. For telescopes within 100 m of 
shower axis it is possible to reject about 80% of the proton 
showers retaining about 70% of 7— ray showers based on flat- 
ness parameter alone. 



5 Conclusions 

In this work we have demonstrated the use of parameters 
based on Cerenkov photon density fluctuations for gamma- 
hadron separation. Using local density fluctuations it is pos- 
sible to reject about 50% of proton showers retaining about 
85% of 7— ray initiated showers. Whereas, based on medium 
range density fluctuations, it is possible to reject about 60- 
70% of proton initiated showers retaining about 80% of show- 
ers produced by 7— rays. Flatness parameter, on the other 
hand, serves as a useful discriminant at core distances away 
from hump. Using this parameter it is possible to reject about 
80% of proton showers, retaining about 70% of 7— ray in- 
duced showers, for core distances within 100 m. Using these 



Table 4. Quality factors with density parameters applied in tandem, 
for core distance < 100 m 



Obs 


Threshold 


Quality factor 


Fraction 


Fraction 


altitude 


value of 




of 


of 


(km) 


LDF, MDF 




accepted 


accepted 




&a 




7— rays 


protons 





0.23, 0.10 


2.221 ± 0.167 


0.655 


0.087 




&0.33 








1 


0.19, 0.09 


2.177 ± 0.175 


0.580 


0.071 




&0.52 








2.2 


0.15, 0.09 


3.365 ± 0.368 


0.563 


0.028 




0.66 
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Fig. 5. Variation of a parameter as a function of core distance for 
showers generated by 500 GeV 7— rays (diamonds and continuous 
line) and 1 TeV protons (+ and dotted line), for altitude of (a) km, 
(b) 1 km and (c) 2.2 km. a parameter is calculated for each sector. 



three parameters in tandem it is possible to improve rejection 

efficiencies further. Table 4 lists the quality factors at vari- 
ous altitudes obtained by applying all the three parameters, 
LDF, MP and a together. Here we restrict to only the tele- 
scopes within the distance of 100 m from shower core. It 
can be seen that it is possible to reject about 95% of proton 
showers, retaining about 60% of 7— ray showers at different 
observation altitudes, using density based parameters in tan- 
dem. In addition, using these parameters in tandem with the 
parameters based on timing information such as timing jitter 
and Cerenkov pulse shape parameters will greatly improve 
the sensitivity of the experiments based on wavefront sam- 
pling technique. 
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Fig. 6. Gamma-hadron separation based on a parameter for sectors 
within 100 m from shower axis for showers initiated by 500 GeV 
7— rays and 1 TeV protons, at observation altitude of 1 km. Dashed 
line indicates the ordinate for optimum quaUty factor. Quality fac- 
tors for this case as well as for other observation altitudes are listed 
in Table 3. 



